Inductive acceleration of UHECRs in sheared relativistic jets 
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Abstract 

Relativistic outflows carrying large scale magnetic fields have large inductive potential and may 
accelerate protons to ultra high energies. We discuss a novel scheme of Ultra-High Energy Cosmic 
Ray (UHECR) acceleration due to drifts in magnetized, cylindrically collimatcd, sheared jets of 
powerful active galaxies (with jet luminosity > 10''^ erg s~^). We point out that a positively charged 
particle carried by such a plasma is in an unstable equilibrium if B • V x v < 0, so that kinetic drift 
along the velocity shear would lead to fast, regular energy gain. This can be achieved in an axially 
inhomogeneous jet through gradient drift induced by propagation of inertial Alfven waves along 
the jet. We show that if a seed of pre-accelerated particles with energy below the ankle < 10^^ cV is 
present, these particles can be boosted to energies above 10^^ eV. A key feature of the mechanism 
is that the highest rigidity (ratio of energy to charge) particles are accelerated most efficiently 
implying the dominance of light nuclei for energies above the ankle in our model: from a mixed 
population of pre-accelerated particle the drift mechanism picks up and boosts protons preferably. 
In addition, after a particle traversed large fraction of the available potential, its Larmor radius 
becomes of the order of the jet thickness. In this case, the maximum possible acceleration rate of 
inverse relativistic gyro-frequency is achieved and a particle finally become unconfined and leave the 
jet. The power-law spectrum of the resulting UHE particles flattens with time and asymptotically 
may become oc E^^. The real injection spectrum depends on the distribution of prc-accclcratcd 
particles inside a jet and, in case of contribution from many sources, on the distribution of total 
potential drop. 

We also point out that astrophysical schemes based on DC-type acceleration (by electric fleld 
parallel to or larger than magnetic fleld) cannot have potentials larger than ~ 10^^ Volts and thus 
fell short by many orders of magnitude to produce UHECRs. 



*Electronic address: lyutikov@purdue.edu 
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I. INTRODUCTION 



Acceleration of Ultra High Energy Cosmic Rays (UHECRs), with energies reaching 3 x 
10^° eV, remains one of the main challenges of modern physics and astrophysics [1-5]. 
UHECRs can be either accelerated by " astrophysical Zevatrons" (or bottom- up scenario) or 
a be a by-product of decay of massive exotic particles (top-down scenario). 

The most commonly discussed astrophysical mechanism of acceleration of UHECRs is 
diffuse shock acceleration by relativistic shocks [6, 7] where small scale magnetic field , 
amplified typically to near-equipartition (see §H), are required. In addition, it is necessary 
that scattering rate is close to the maximum (Bohm) rate in order to overcome advection of 
particles downstream of a shock. It is far from clear whether such rates are achievable [8] . 

Commonly, it is believed that low energy cosmic rays , with energies approximately be- 
low the "ankle" at ~ 10^^ eV, are of Galactic origin (GCR), while those above the "ankle" 
are extragalactic. [Note, that exact location of transition between Galactic and extragalac- 
tic cosmic rays is not clear and may instead be associated with the "second knee", [9].] 
Since spectra below and above the "ankle" are different, different acceleration mechanisms 
should operate [see 10, for an alternative holistic scenario]. In this paper we propose new 
mechanism that operates at highest energies and thus is responsible for extragalactic cos- 
mic rays (EGCR): particle acceleration by inductive electric fields in magnetized, relativistic 
sheared outflows. 

Inductive electric fields and cosmic ray acceleration has been mentioned by many authors 
[1, 11-13], but no concrete acceleration mechanism has been proposed (with a possible 
exception of [14], see below). It was only pointed out that there is a large electric potential 
and it was assumed thai particle somehow manage to tap it. Since inductive electric fields are 
orthogonal to magnetic field and particles cannot move freely along them, it is not obvious 
how to achieve energy gain. Another acceleration scheme is acceleration in DC electric field 
aligned with magnetic fields (double layers), but these, as we demonstrate in appendix A, 
cannot account for UHECRs. 
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A. Inductive electric fields 



As mentioned above, inductive electric fields are not easily accessible for acceleration since 
particle need to move across magnetic field , a processes prohibited under ideal Magneto- 
Hydrodynamics (MHD) approximation. On the other hand, kinetic drifts may result in 
regular motion across magnetic field and along electric field leading to regular energy gain 
as compared to the stochastic, Fermi-type schemes. Since the drift velocity increases with 
particle energy, the rate of energy gain will also increases with particle energy. Thus, the 
highest energy particles will be accelerated most efficiently. This means that from a pre- 
accelerated population, the mechanism proposed here will pick up particles with highest 
energy and will boost them to even higher values. In addition, when a particle has crossed 
a considerable fraction of the total available potential, it gyro-radius becomes comparable 
with flow scale (§11). In this case, drift approximation brakes down. As a result, as long as 
the particle remains inside a jet, the acceleration rate reaches the theoretical maximum of 
an inverse gyro-frequency (^VB). Thus, the most efficient acceleration occurs right before 
the particle leaves the jet. One may say that in case of inductive acceleration, becoming 
unbound is beneficial to acceleration, contrary to the case of stochastic acceleration when for 
unbound particles acceleration ceases. 

Typically, the direction of drift is along the normal to the magnetic field and to the di- 
rection of the force that induces a drift. In sheared cylindrical jet with toroidal magnetic 
field electric field is in radial direction, so that in order to gain energy particle should ex- 
perience radial drift. It is then required that there should be a force along the axis. Such 
force may arise if a jet is axially inhomogeneous resulting in gradient drift due to changing 
magnetic pressure. Several types of inhomogeneity can occur in astrophysical jets. Nu- 
merical simulations of hydromagnetic jets show complicated axial structure with regions of 
compression and rarefaction [15-17]. Besides shocks, these simulations show prominent sub- 
sonic density and magnetic field variations inside jets. Subsonic variations propagate roughly 
with the jet fiow and induce axial variations of magnetic field leading to radial drifts. An 
important property of the mechanism proposed here is that for a particle at rest with respect 
to the bulk flow energy gain (or loss) is independent of the direction of the drift {e.g. away 
or towards the axis; see §V). On the other hand, if an axial variation is periodic, e.g. a wave, 
changing the direction of the drift would lead to alternative energy gains and losses. If in 
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one wave period a high energy particle acquires a large enough energy so that its gyro-radius 
becomes comparable to jet scale, it will leave the jet. 

B. Electric drifts and particle acceleration 

To separate effects of electric drift in crossed electric and magnetic fields (which does not 
lead to energy gain) and particle acceleration, we should consider particle dynamics in a local 
plasma rest frame, where electric field vanishes at the origin. Since the velocity field is likely 
to be sheared {e.g. higher velocity near the axis of a jet), there will be non- vanishing electric 
field away from the origin. Value of the electric field at a position of a particle will be smaller 
than inductive field in the laboratory frame, approximately, by a ratio of Larmor radius to 
shear scale. Typically, both shear velocity and electric field would be linear function near 
the origin, so that electric potential is a quadratic function of distance from the origin of the 
local plasma rest frame. Depending on the sign of the coefficient in front of the quadratic 
term, electric potential has either a minimum or maximum near the origin. This depends 
on the sign of the charge and sign of the product B • V x v. In §111 we show that in sheared 
flows with B • V X V < protons are near the maximum of electric potential, so that drift 
motion away from the origin {along or against the direction of the shear), will lead to energy 
gain. 

Our model bears some resemblance to Bell model [14] of cosmic ray acceleration at pulsar 
wind nebular, which relies on large (unreasonably large, in our view) magnetic field gradients 
near the axis of the expanding fiow in order to "pull-in" pre- accelerated charges against 
weakly relativistic expansion of plerion. The proposed mechanisms also has both similarities 
and differences from the shock drift acceleration at non-parallel shocks [18]. A shock can be 
considered as an extreme case of field inhomogeneity (a delta-function), so that drift velocity 
is a delta function [19]. This can similarly produce energy gains due to displacements of 
particle along inductive electric field associated with shock motion. The idea that particle 
energy can be hmited by "lateral" escape, is also reminiscent of Eichler model [20] of particle 
acceleration at the Earth's bow shock. The difference is that in Eichler's model the sideways 
motion of a CR was due to cross-field diffusion, while in our case it is a gradient drift. 

Our acceleration scheme is also very different from those previously considered in sheared 
fiows due to second-order Fermi process [3, 21-23]. In the references cited above particle 
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energetization occurs due elastic scattering by magnetic field inhomogeneities in the flow 
rest frame. Acceleration is stochastic, while acceleration time scales are typically very long, 
scaling as velocity shear squared, V^, [22], so that extremely narrow layers are required. 
Contrary to this, in our scheme acceleration is regular and scales as V. 

Let us also point out that any acceleration scheme based on DC-type acceleration, where 
there is large electric field along magnetic field or where electric field becomes larger than 
magnetic field, cannot accelerate particles to energies above 10^^ eV, appendix A. This 
excludes, for example, magnetic reconnection [24, 25]. 

In this paper we start in Section 2 by briefly describing on general basis the properties 
of UHECR acceleration sites. The constraints on source energetics and size arc discussed 
within the context of cylindrical AGN jets. In Section 3 and 4 we study the evolution of 
particles trapped in sheared flows and in the presence of Alfven waves. The details of the 
acceleration mechanisms at play are described in details in Section 5. The corresponding 
particle spectrum is derived in Section 6 before considering the astrophysical implications 
in Section 7. A discussion is presented in Section 8 and a brief conclusion in Section 9. 

II. AGN JETS AS POSSIBLE SITES OF UHECR ACCELERATION 
A. Luminosity - Potential relation 

In bottom-up acceleration schemes [4, 26] any source that accelerates a particle to a 
particular energy should have a total electric potential drop of the same order (it is unlikely 
that a particle completes many circuits with much smaller potential drops). On very basic 
grounds, we expect that in any astrophysical system electric field are no larger than magnetic 
field . [Locahzed gaps, regions with either parallel electric field or total electric field larger 
than magnetic, may exist but the typical potential drop in them is small (for purposes 
of UHECR acceleration) and is limited by lepto-photonic pair production, [27].] This has 
two consequences. First, regardless of the acceleration mechanism, the highest possible 
acceleration rate is inverse of the relativistic gyro-frequency, ~ ^/Zcub where cub — eB/mc 
is cyclotron frequency, Z is a charge of a particle in terms of elementary charge e, m is a 
mass of a particle. Secondly, we can place very general constraints on possible accelerating 
sites [see also 28]. If we parameterize E = j3oB where E is the electric field and B is the 
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magnetic field (/3o is the ratio of typical electric to magnetic fields ) and if an accelerating 
region has a size R the available potential is $ ~ RE ~ PoRB and the energy of particle 
crossing that potential would be 

S = ZePoBR (1) 

For /5o ~ 1 this is also the condition that the particle's gyro-radius is smaller than R [29]. 
Thus, for /3o ~ 1, a system will be able to accelerate particles to the maximum energy 
that it can confine. We would like to stress that the conditions on acceleration of UHECRs 
(1) are different from confinement condition ri < R, where ri is the Larmor radius. The 
two conditions are different by a factor /Sq: for f3o 1 a system can confine particles with 
energies much larger than those to which it can accelerate. 

Since BR is a measure of the total current in the system, Eq. (1) can be rewritten as 
(assuming /3o ~ 1) 

$ = /3of (2) 

This determines the maximum electric potential and corresponding energy S — Ze^ to 
which a system carrying a current / can accelerate. 

Unless the electric field is strictly along the magnetic field (e. 5'. in double layers, see 
appendix A), there is Poynting flux in the system. Since Poynting flux cannot be larger 
than the total luminosity of a given source, assuming a spherical expansion we find 

L = 4rf — c R^B^f3oc = — 3 

a Att a cr po 

where a is the ratio of Poynting to particle flux 

Thus, assuming /5o ~ 1, 



L = i±^/* (4) 

a 

where we introduced coefficients of the order of unity so that in the limit cr — 00 the system 
impedance is that of free space, Att/c; L^q — erg/s and for numerical estimate we 

assumed a — 1. For each particular system these relations will be modified by a numeri- 
cal factor which depends on geometric properties of the system. The case of cyhndrically 
coUimated AGN jets is considered in §11 B. 



Eq. (4) relates the maximum available potential to the luminosity of the source. Assuming 
that UHECRs are protons, we can select possible candidates for UHECRs acceleration based 
on their luminosity: to achieve 3 x 10^*^ eV we need L > 10^^ erg/s. We can immediately 
exclude acceleration of UHECRs in low power AGNs (e.g. Cen A and M87), low power BL 
Lacs and starhurst galaxies (e.g. M82 & NGC 253) (in case of proton component), see also 
§VII. This limits the possibilities to more distant high power AGNs like higher power FR I, 
FR II radiogalaxies, radio-loud quasars and GRBs. [In appendix B we argue against GRB 
shocks as source of UHECRs.] Luminosity constraint (4) also excludes, in our view, hohstic 
scenario, [10], where the same acceleration mechanism operates in all galaxies, including our 
own. Only "special'' galaxies, satisfying condition (4) are able to accelerate UHECRs. Note, 
that in order to maximize potential at a given luminosity, it is required that cr > 1, but not 
necessarily ^ 1. 

Another constraint that acceleration sites should satisfy is that radiative losses should 
not degrade particle energy. We can derive very general constraints on possible location of 
cosmic ray acceleration just by balancing most efficient acceleration, by E r-j B, and radiative 
losses [see also 31]. Typically, radiative losses scale as oc j^u where 7 is the Lorentz factor 
of a cosmic ray and u is energy density of the medium that leads to radiative losses. For 
synchrotron losses u — ub — B'^/8tt, for IC losses u — uic is the energy density of soft 
photons. If we normalize total energy density to energy density of magnetic field = C'^Bi 
C > 1, and equate the maximal rate of energy gain d£/dt — ZecE ~ ZecB to radiative 
losses dS/dt ~ ^^u)\{E j mc^Y we find 



where Tc = jmc?. This may be considered as a second constraint on the site of particle 
acceleration. Resolving eqns (1) and (5) with respect to B and R we find 



Where we also allowed a possibility that plasma is expanding with bulk Lorentz factor F, so 
that magnetic field in the plasma rest- frame is BjV and typical scale is i?/F, and assumed 




(5) 





C = 1; EeV = 10^^ eV. 
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Relations (6) show that higher energy cosmic rays are better accelerated at large dis- 
tances. AGN jets, which propagate to more than 100 kpc distances present an interesting 
possibility. Note, that as long as the jet remains relativistic, the total inductive potential 
is approximately conserved, so one can "wait" a long time for a particle to get accelerated 
without worrying about radiative loses. Thus, UHECRs can be accelerated inside the jet at 
distances from a fraction of a parsec (Eq. (6)) to hundreds of kpc, as long as the jet remains 
relativistic and sustains a large inductive potential. 



B. Cylindrical AGN jets 

As an illustration of the general principles discussed above let us consider force-free 
cylindrically colhmated jets. A possible structure of such a jet is a relativistic sheared 
diffuse pinch 

1 + {r/RcY 
Er = ^1 - l/r{rfB^ 

^ l + {r/R^r^' 

where Rc is a radius of a current-carrying core. System remains in force-balance for arbitrary 
r(r) [32]. Alternatively, instead of Bz particle pressure may ensure force-balance. 
For illustration purposes let us choose a Lorentz factor with a profile. 



r(r) = 



where Fc — 1/ ^Jl — (51 is the Lorentz factor of the fast core and we assumed that Lorentz 
factor changes on the scale of the core radius. Then the total potential is 

BR 2 

$tot = -Boi?carctanh(/3e) ~ — ^ In , (9) 

while total the Poynting flux is 

^ = — (-a^-<=t-i'(A) - ft j - — T3ft (10) 

where the last equahties assume /3c — > 1- Thus 

*« = /f\^^-/fV^ (11) 



If a particle start on the axis with 7 = 7c, and drifts in a radial direction, its Lorentz 
factor 7 changes with radius as 

7 = 7^+ ^^'"^"^ |^arctanh/3c - arctanh ^ ^ ^^^"^ ^ j (12) 

so that the particle gains large fraction of the potential after reaching r ~ RcTc- At this 
point its Larmor radius would become 

~ arctanh/3c > 1 (13) 

Rc^c 

This illustrates an important point: after a particle has crossed a large fraction of the total 
available potential, its Larmor radius becomes comparable to the flow scale, so that it will 
decouple from the flow. 



III. POTENTIAL ENERGY OF A CHARGE IN SHEARED FLOW 

Our model of UHECR acceleration relies on the observation that in a transversely sheared 
flow one sign of charges is located at a maximum of electric potential, as we describe in this 
section. Consider sheared flow carrying magnetic flcld . At each point there is electric field 
E = —V X B/c, so that the electric potential is determined by 

A$ = -V • (v X B) = - (B • (V X v) - V • (V X B)) (14) 
c c 

In a local rest frame, the second term in Eq. (14) vanishes at the position of a particle and 
is generally sub-dominant to the first term for v <^ c [84]. Then V x B = and we find [85] 

A$ = - (B • V X v) (15) 
c 

We have arrived at an important result: depending on the sign of the quantity (B • V x v) 
(which is a scalar) charges of one sign are near potential minimum, while those with the 
opposite sign are near potential maximum. Since electric field is perpendicular both to 
velocity and magnetic field, locally, the electric potential is a function of only one coordinate 
along this direction. For (B • V x v) < (we will call this case negative shear) ions are near 
potential maximum. 

The procedure outlined above to calculate electric potential is beyond the limits of ap- 
plicability of non-relativistic MHD, which assumes quasi-neutrality and thus neglects the 
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dynamical effects associated witfi tfie potential (15). Thus, even in the low frequency regime 
with non-relativistic velocities, a conventional realm of MHD theory, one should use at least 
two fluid approach and also must retain both charge density as well as displacement current 

in Maxwell equations. 

We can arrive at a more general result using a relativistic MHD formulation. An invariant 
charge density is given by 

(16) 

where j^^ is four-current, is fluid four-velocity and F^^" is the electromagnetic tensor and 
Greek indexes run 0,1,2,3 and we set speed of light to unity here. In the framework of 
relativistic MHD, when stress-energy tensor is diagonalizable, the electromagnetic tensor 
can be written as F^^'^ — e'^^'^^BaUp, where Ba is a four-vector of magnetic field and e^^°'^ is 
a unit antisymmetric tensor. Charge density becomes 

P = -^e'^'^"^^^ {B^,pu(, + B^up,,) (17) 
47r 

If the magnetic field is stationary and homogeneous, Ba,f3 = 0, then we arrive at 

P=^^'^X>- (18) 

which generalizes Eq. (15) to an invariant form. 

Under ideal fluid approximation particles cannot move across magnetic field lines, so that 
they cannot "sample" the electric potential (15). On the other hand, kinetic effects, like 
drift motions, may lead to regular radial displacement along the shear and thus along electric 
field. In this case one sign of charge will be gaining energy, while another sign will be losing 
energy. This is independent on whether the drift is along the shear or counter to the shear 
direction and thus is independent on the sign of the magnetic field gradient that induces the 
shear. 

For low energy particles, with small but finite Larmor radii, variations of potential energy 
need to be taken into account for a correct description of the plasma distribution function 
in magnetized flows with transverse shear [34-36]. [Since the distribution function should 
depend on constants of motion, in presence of a velocity shear these constants are different 
from the case of no shear.] Kinetic motion allows particle to probe regions with different 
potential and lead to redistribution of positive and negative changes in an effort to screen 
the shear-induced electric field. These effects are absent if background plasma is cold. On 
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the other hand, for a small population of high energy particles with large Larmor radii, 
electrostatic potential (15) may become important. 

IV. ACCELERATION MECHANISM IN A NUTSHELL 

Our key point is that an electrostatic cross-field potential (15) may be used for particle 
acceleration. In this Section we briefly outline the principal astrophysical and plasma phys- 
ical issues related to this acceleration mechanism. In appendix C and §V we elaborate in 
more details the qualitative statements made in this Section. 

There is a consensus that relativistic outflows, and AGN jets in particular, are acceler- 
ated to relativistic speeds and coUimated by large scale magnetic fields threading accretion 
disk and central black hole. This is based mostly on theoretical and numerical models of 
jet launching from black hole- accretion disk systems [37-40]. In addition, polarization 
properties of parsec scale AGNs jets [32] and, possibly, GRBs [41] are also consistent with 
the presence of large scale magnetic field . Energetically, magnetic fields may carry a large 
fraction of jet luminosity [42, 43], so that the parameter a introduced above is not too small, 
(J > 1. 

As the jet expands, the ratio of toroidal to poloidal magnetic field increases approximately 
linearly with cylindrical radius of a jet, so that at largest scales magnetic field is dominated 
by a toroidal component B^. In addition, toroidal magnetic field may provide jet coUimation, 
so that asymptotically AGN jets may also be fully coUimated to a cyhndrical shape [44]. 
One also expects that jets are sheared, so that the central spine of the jet is moving with 
larger velocity than its periphery [45-47], so that = Vz{r) (r is the radial coordinate in 
cylindrical geometry) is decreasing outward. In this case, (B ■ V x v) = —B^d-r-Vz- If, in 
addition, the jet is inhomogcncous along the axis, B^p = B^[z), there will be a gradient drift 
along r, which in case of negative shear would lead to proton acceleration. In another type 
of jets, with positive shear (B • V x v) > 0, electrons are accelerated while protons lose their 
gyrational energy eventually coming to a rest with respect to the bulk flow. 

An important requirement on the proposed acceleration scheme is that the flow be nearly 
cylindrically coUimated. This stems from the fact that as long as a particle has not crossed 
the total available potential, its drift velocity is non-relativistic, so that in conically expand- 
ing relativistic flows energy gain cannot compete with adiabatic losses. 
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Next we discuss the salient properties of kinetic particle drift motion and energy gain 
in sheared flows. Consider a plane-parallel sheared flow in which the magnetic field is per- 
pendicular to the velocity. For definiteness, let magnetic field be directed along x axis, 
= Bq and velocity along z axis, V = Vz{y) (see Fig. 2). Let us transform to 
a frame moving with the fluid at some particular values of y = yo, e.g. at yo = 0. 
Expanding velocity near y — 0, Vz — Pocy/Ly, where is characteristic velocity in 
terms of speed of hght and Ly is a characteristic scale. Thus, there is electric held 
E = —VzBx/cCy — —{Vo/c)Boy/Lvey — —PoBoy/LvCy and electric potential 

^ = PoBo:^ (19) 

Next, assume that there is a long wavelength inertial Alfven wave propagating along 
the 2;-direction with a phase speed Va- For Va <^ c, the magnetic perturbation SB in the 
wave is much larger than electric perturbations 5E by a factor c/Va, so that the wave is 
nearly magneto-static (in other words cu <S kzc). A test particle will experience a drift in 
the x-direction with magnitude 

6B -fc'^k, -fc^k, 
Ud ~ —- (20) 

where we assumed strong perturbation 5Bq B and ujb,q — eBo/mc. As a particle drifts 
along electric fleld in the x-direction its Lorentz factor evolves according to 

,2+ „,2^3k24 



ZeEyUd ZujBfiuit Yc^kjt 
mc^ cLy Zujb o^v 



Thus, energy gain or loss of a test particle depends on direction of magnetic field , sign of 
charge (through cub) and direction of velocity vorticity (through sign of Ly): in other words 
it depends on sign of shear. On the other hand, it is independent of the direction of the 
drift. 

Since drift velocity increases with particle energy, the rate of energy gain also increases 
with energy, see Eq. (21). This leads to one of the most unusual properties of the proposed 
acceleration mechanisms: highest energy (or highest rigidity) particles are accelerated most 
efficiently. In addition, at the last stages of acceleration, when particle Larmor radius 
becomes of the order of jet scale, particle motion in positive shear fiow becomes unstable 
even without gradient drift while acceleration rate does reach absolute theoretical maximum 
of inverses relativistic gyro-frequency (§VB). 
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In addition, since acceleration rate is inversely proportional to charge, at a given energy 
small charge (higher rigidity) particles are accelerated most efficiently. Thus, from a pop- 
ulation of pre-accelerated particles with mixed composition, drift mechanism will pick up 
particles with smallest charge: protons. This explains why above the ankle protons start to 
dominate over heavy nuclei. 

Above is the upshot of the proposed acceleration mechanism. Next we elaborate its de- 
tails. In appendix C we consider properties of electromagnetic waves in transversely sheared 
cold plasma; this is required to in order to set up correctly the background electromagnetic 
fields for the drift motion of a particle. In §V we consider drift motion of test particle in 
an electromagnetic fields of a low frequency inertial Alfven wave propagating in a sheared 
fiow in case of weak shear, when particle gyro-radius is much smaller than shear scale. This 
is done analytically for weak waves and numerically for strong waves. In §VB we consider 
particle motion for strong shear, when shear scale is comparable to gyro-scale, so that drift 
approximation is inapplicable. In this case a new type of acceleration mechanism turns on, 
which does not require presence of velocity-aligned magnetic field gradients and which leads 
to theoretically maximal acceleration rate on inverse of relativistic gyro-frequency. 

V. DETAILS OF THE ACCELERATION MECHANISM 
A. Stage 1: Kinetic drift 

1. Kinetic drift in linear Alfven waves 

In this Section we first consider kinetic motion of a particle in a combined electromagnetic 
fields of sheared fiow and a low frequency inertial Alfven wave propagating along velocity. 
For simplicity we concentrate on planar problem and assume a linear velocity profile. Using 
results of appendix C, electromagnetic fields of the inertial Alfven are then given by 



SEy = aBo sin {k,{z - {Va - V,)t)) 




(22) 
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where a is a non-linearity parameter of the wave, defined here as a ratio of electric field 
amplitude to background magnetic field . Assuming linear shear, Vz — rjy, Eq. (22) gives 

Ey = Bo (-^ + asmkz{z - {Va - w)^) 
B, = Bo(l- sin kz{z - {Va - vy)t)] (23) 

where rj = V and we assumed that a test particle is located near the zero wave phase. 
The equations of motion of a test particle are 

dt {^d,v) = ^ (e + ^ X (24) 
Next, we separate particle motion into slow guiding center motion R(t) and fast rotation 

m 

r = Ro{t)+m (25) 

with 1^1 of the order of Larmor radius ~ r^. Assuming long wavelength perturbations A;^ ^ 0, 
weak perturbing wave a ^ 0, and expanding in small gyro radii ri = \J'yo — 1/ujb (we 
neglect motion along magnetic field, = 0), we find equations of motion of the gyrations 
center. 



dtYo = -^dtZo + - k,t{Zu:BViho - v) 



70 70 2 ^ V - - "2BoVa 



2^ f Zujb \ c... . , SE 



where capital letters denote coordinates of Rq. Equations (26) can be integrated to give a 
drift in the y-direction with velocity 

Ud _ c fc,c7o 1 - (Va/c:)^ + V1o/{Zu;b) c(1 - {VA/cf) k.cjo 

— (X — r ^ d V J 

c Va2Zub 1 - Tj'yo/iZujB) Va 2Zub 

where we have reinstated the speed of light explicitly and the last approximation assumes 
small 77 — > 0. In addition, there is insignificant motion along z which results in a 2; displace- 
ment of the order of k^Sz ~ a{c/VA) < 1- 

The origin of terms in Eq. (27) can be easily understood: a{c/VA) is a relative amplitude 
of magnetic field perturbations due to presence of the wave, kzC'-)Q/2u}B is a ratio of particle 
Larmor radius to inhomogeneity scale and the last terms represent relativistic and shear 
corrections. Note, that as plasma becomes strongly magnetized, Va ^ c, value of the drift 
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velocity decreases. This can be understood by noting that for subluminal waves there is 
always a frame where electric field is zero and magnetic field is constant [48]. What is 
important for drift motion is the value of the magnetic field in that frame. For shear Alfven 
waves this frame is moving with the velocity V^. For Va c, the value of magnetic field in 
this frame decreases which leads to smaller transverse drift velocity. 

For order-of-magnitude estimates, we can assume Va ~ c. [Note that for -C c, for 
a given potential total energy fiux increases by (Va/c)^, while for Va — > c drift velocity 
becomes small, so that the mechanism operates best at Va < c] Then 

- - (28) 
c 2Zu;b 

The energy of a particle changes both due to shear-induced electric field as well as the 
electric field of the perturbing wave. We are interested in the former contribution; the 
latter process - gradient drift along the electric field of the wave gives a well known surfing 
mechanism of particle acceleration [49]. Energy gain in such a process is limited by electric 
field of the wave and spatial dimension of the wave. Thus, if wave amplitude is of the order of 
unity, 6B/B ~ 1 wave surfing will lead to similar energy gains as drift in a sheared flow. In 
astrophysical applications we expect 6B/B < 1, so that shear-induced potential is typically 
larger than the electric potential associated with perturbing wave; see also §V A 2 for more 
detailed discussion of effects of surfing. 

Assuming that a particle keeps coherence with the wave (see §VA2), its energy in a 
shear-induced electric field evolves according to 

where we introduced inhomogeneity scale of magnetic field Lb = i/kz and shear scale 
Ly ~ c/V — c/t). Note that variations of energy are independent of the sign of magnetic 
field gradient Lb, but depend on particle charge, magnetic field and electric field gradient. 
Eq. (29) implies very fast energy growth: 

where 



2v^Lb / Lv 
To = \ —. — (31 
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Note that the acceleration time scale is inversely proportional to the square root of initial 
particle energy: 

VTo 



Thus, more energetic particles with smaller charge are accelerated on shorter time scales. 

Initially the energy grows quadratically with time. At later times, the Larmor radius of a 
particle starts to increase in response to increasing energy. This leads to faster drift and faster 
rate of energy gain, which formally becomes super-exponential (finite time singularity). Of 
course, in order to get infinite energy a particle should move infinitely far, so that in reality 
maximum energy gain will be limited by the size of the system and it's total available 
potential. In addition, before the finite time singularity is reached, drift approximation 
will break down. This will occur when the particle Larmor radius becomes of the order of 
inhomogeneity scales L^, Le- This happens only a short time before finite time singularity, 
at times 



"Tacc {^^) 



2. Nonlinear waves: effects of wave surfing on energy gain 

We have shown that if a particle can keep in phase with the wave in a combined elec- 
tromagnetic fields of sheared fiow and inertial Alfven wave, its energy increases at super- 
exponential rate. In case of a linear wave this assumption is not completely satisfied: the 



acceleration time scale (32) is longer than the wave period by a factor yLy/r/^. This dif- 
ficulty may be overcome if perturbing wave is non-linear, with finite relative amplitude 
a — SE/Bq. In this case the particle will also experience a drift in crossed and Bq fields. 
If a particle is located at such phase of the wave so that this drift velocity is along the di- 
rection of propagation, it may remain in resonance for a long time. For substantially strong 
waves the longitudinal drift velocity (along the direction of wave propagation) ~ 5E/Bq 
may become of the order of the phase velocity of the waves. Such particles arc effectively 
trapped by the wave. In surfing acceleration schemes this leads to particle energy gain as 
it drifts along the electric field of the wave. In our case trapping leads to nearly constant 
velocity drift along the shear-induced electric field . 

To test this possibility we consider numerically particles motion in a sheared fiow with 
finite amplitude inertial Alfven wave. As a test case, we consider planar problem, similar 
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to the one described in Fig. 2. For background velocity profile we chose — VAta,nh{y/ys), 
where is a parameter determining the strength of the shear. On top of the shear we impose 
an inertial Alfven wave, with fluctuating flelds given in the WKB approximation (appendix 
C) by Eqns (22). We then seed the wave with test particle at random phases in the wave 
and follow their trajectory. Results of simulations are presented in Fig. 5 where we plot the 
maximal energy gained by a particle in a non-linear inertial Alfven wave propagating in a 
sheared flow as a function on initial phase for the following parameters: B — B^, Va — 0.5c, 
a — 0.25, kz — 0.025, = —5 (negative is needed to have a negative shear; for positive 
there is no acceleration), initial energy 70 = v^, initial momentum is directed along y axis. 
Upper curve corresponds to the case of shear shear flow and an inertial Alfven wave, lower 
curve is the case without shear. Fig. 5 shows that particle indeed can gain large energies in 
a sheared flow: most of the energy gain is due to motion in sheared electric field and not due 
to wave surfing: the wave mostly provides a drift velocity along the shear. Without a wave, 
particle just executes gyrational motion near the point y = 0. 

Another related type of acceleration in strong waves should be mentioned. In a strong 
inertial Alfven wave propagating across magnetic held , fluctuating magnetic field of the 
wave 5Bx at some wave phases reduces total magnetic field , so that total electric field 
(which is a sum of wave electric field and drift-induced field) may become larger than the 
total magnetic field . In this case, there exists a reference frame where the magnetic field is 
zero and only electric field is present, so that particles will experience strong acceleration. 
This type of acceleration is similar to surfing acceleration in strong electrostatic (upper 
hybrid) waves [49]. Since both charges of particles are accelerated in opposite direction, 
such mechanism cannot be responsible for UHECRs, see appendix A. For the parameters 
of our simulations condition E < B is always satisfied. 

B. Stage 2: large Larmor radius 

When the Larmor radius becomes comparable to shear scale particle dynamics changes 
qualitatively. As we show below [see also 34], in case of narrow shear layers, when Lar- 
mor radius becomes comparable to shear scale, particle motion becomes unstable even for 
homogeneous flow, without perturbing inertial Alfven wave. 

The easiest way to consider motion of a test particle in such a case is using Hamiltonian 
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formalism. Let's assume that magnetic field is in the x-direction and that there is a shear 
velocity along z with v^iy) — P{y)c (Fig. 2) and transform to the frame where v{0) — 0. 



Hamiltonian of a particle in a sheared flow, if = y {mc^Y + c^{P — eA/c)^ + e$, where P 
is generalized momentum, A is vector potential and $ if electric potential, can be written 
as 

^_4J,^4^ + ?M(pZ^£^) (34) 

where Py is momentum along y and Yg = Pz/muoB is y coordinate of the guiding center. For 
linear shear the Hamiltonian is 



2c2 



with the corresponding equations of motion 

Py/m 



dty = 



V m2c2 ^ c2 



{Yi, — tj)Z'^uJi,m 

OtPy = -yrjZujBm H — , ' ' = (36) 



Let us assume that center of gyration is at Yg — 0. After dimensionalizing, py/mc — > Py, 
y — > cy/ZujBi t — > cuBt, rj — > rjcuB, Eq. (36) becomes 

Py 



dty = 



^Jl+Pl + y' 



From conservation of energy 70 = y 1 + + y^ + riy /2, we find that the turn-around point, 
where p^^ = is at ymax — V^^/l + 7o^ + a/I + 27o?7 + rj'^/r). In addition, there exists a 
critical trajectory, so that as time goes to infinity a particle moves with constant velocity 
along z axis. In this case dtPy — and we find ymax — y'l/jy^ — 1. Equating the two 
expressions for ymax we find critical values of shear 

Vcrit = Vto - 1 - 7o ~ -77- for 70 > 1 (38) 

^7o 

and ymax,crit = Y^2/3o/(l - /3o), where /3o = ^1 - I/70. Thus, for 77 < r]crit < 0, orbits are 
unbound and energy growth exponentially. In the limit 70 ^ 1 this occurs when the shear 
scale Ly — c/V is half of gyro-radius, Ly — c^q/2ijJb- 
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Particle trajectories can be found in quadratures: 



4(l + y2) 
\ (27o - y^^)^ 



t 

(see Fig. (3)). In the non-relativistic limit equations of motion can be integrated exactly, 



y = tlcos Z(jJb\Ji + rj/ZujBt., z = , ^ — sin Zcu^^/l + r^/cugt (40) 

^Jl + r]/Zu;B 

This clearly shows that for strong negative shear, rj < —ujb, particle trajectory is unstable 
and its energy growth exponentially. For positive shear, rj > 0, particle motion is stable. 
Thus, when measured in term of particle gyration frequency cu/^q, critical shear reduces 
from —1 in the non-relativistic limit to —1/2 in the ultra- relativistic limit. 

Similar results may be obtained numerically for velocity profile oc tanh(y/ys). This 
case has an advantage over the linear profile since condition E < B is satisfied everywhere. 
For any given 70 there is a critical value of ys,crit < 0, so that for ys,crit < Vs < ^ particle 
trajectories are unstable, see Fig. 4. 

Energy gain occurs on time scale of Ly/cjo, which for strong shear is of the order of rela- 
tivistic Larmor frequency. Thus, in this case acceleration rate reaches the absolute theoretical 
maximum. 

This acceleration mechanism can be contrasted with shock acceleration, in which case 
acceleration stops when particle Larmor radius becomes of the order of system size. In 
contrast, inductive acceleration become more efficient when particle is unconfined. 



VI. SPECTRUM OF ACCELERATED PARTICLES 

The acceleration mechanism proposed here requires a presence of pre-accelerated seed of 
particle with large Larmor radius. Pre-acceleration mechanism may operate outside of the 
jet; radial drift then may "pull" particles inside a jet. Since magnetic field at the periphery is 

relatively small (e.g. falling off as 1/r), even relatively small energy particle would have large 
Larmor circle and will experience fast radial drift. This can serve as an injection mechanism 
if at the base of the jet, near the central galaxy, there is a population of pre-accelerated 
particles. 

In a steady state, the total particle energy gain is determined by the local electric poten- 
tial, which, in turn, depends on magnetic field and velocity distribution. These are highly 
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uncertain functions for the AGN jets. As an alternative model problem we consider tem- 
poral evolution of distribution of pre-accelerated particles injected into jet at initial time 
and evolving due to gradient drift in sheared flow. Assuming 7^1, the evolution of the 
distribution function (/ = ^) in the rest frame is described by 

dtf + d,{dt{^)f)=S (41) 

where S is density of sources. Assuming impulsive injection of seed population, integration 
along characteristics (30) gives a solution 

/(7,t) = /o(7o(7,t),t = to)(^) (42) 
If the initial injection spectrum is power- law, /o oc I/70, then 

/(7,t)ocl(l+,(l)y" (43) 

~ 7"^ for 7tV'^o > 1 

Thus, for p > 2 the spectrum flattens with time. 

The hardest spectrum that can be achieved has a power law index of 2. This limiting 
case corresponds to unlimited acceleration in a plane-parallel geometry, which is realistically 
applicable to energies well below the total available potential. At highest energies the final 
spectrum will depend on the distribution of pre-accelerated particles with respect to the 
electric potential and, in case of contribution from many sources, on distribution of total 
potentials. 

Relation between the observed spectrum and spectrum at the source is not simple and 
depends on source distribution and (poorly constrained) extragalactic magnetic fields . Since 
high power AGNs are located at the edge of the GZK sphere, a flat spectrum at the source, 
with p approaching 2, seems to be required in order to account for AG AS A observations 
of UHECRs above the ankle [2, 50]. We note however that in the high energy region we 
discuss here currently there seems to be a disagreement between the AGASA ground array 
([51]) and the HiRes fluorescence detector ([52]). Much larger experiments such as Auger, 
Extreme Universe Space Observatory (EUSO), and Orbiting Wide-angle Light- Collectors 
(OWL) will put much better constraints on the injection spectrum. 
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VII. ASTROPHYSICAL VIABILITY 



One of the key observational test of UHECRs origin is the identification of the transition 
from the galactic to extragalactic components. The two competing models are the "ankle 
scenario" [see, e.g. , 53, for a general discussion of the transition], when the transition is 
identified with the intersection of a steep galactic and flat extragalactic spectra, and the 
"dip" scenario, when the transition is identified with the "second knee" [see, e.g. , 54]. Our 
model is consistent with the dip scenario as we discuss below. 

Chemical composition In our model the EGCRs and the GCRs have different spatial 
and astrophysical origin. We assume however that the mechanism driving the pre- accelerated 
component at the acceleration cite of UHECRs is also at play in our galaxy. Prom a pool of 
pre-accelerated particles, the ones with the highest rigidity are accelerated most efficiently. 
Thus, the population of pre-accelerated CRs is a mixture of heavy and light nuclei, the drift 
mechanism will preferably accelerate protons. Our model thus predicts an increase in proton 
over heavy elements as the energy of the UHECR increases. Such trends have apparently 
been observed in the composition of UHECRs with energies between ~ 5 x 10^^ eV ("second 
knee") and 10^^ eV and extending to energies above 10^^-^ eV ([65, 66]). [See [67] for a 
discussion of systematic differences between different measurements and how it limits our 
knowledge of the composition of CRs of the highest energies.] 

Injection spectrum The hardest injection spectrum the model can account for is power 
law index of 2. The dip models generally prefer softer spectra, with index 2.5—2.7; such softer 
spectrum can be achieved both due to spread out distribution of pre-accelerated particles 
with respect to the total potential inside a given jet and by a combination of contribution 
from many sources with different potential drops [54, 68]. For the case of injection spectrum 
with index 2 the distribution of the sources with respect to maximum acceleration energy 
Emax required to explains the observational data is dn/dEmax oCmax [68]. The model 

also naturally provides a lower cut-off to the injection spectrum (see below). 

Density of sources: isotropy and clustering of UHECRs. We associate accel- 
eration cites of UHECRs with powerful FRII galaxies. Assuming proton composition, to 
account for energies above 10^*^ eV the source luminosity should be no less than 10^^ ergs/s 
(sec §2). Assuming a ~ 1, this translates into a limit on a total power of a jet. The total jet 
power, as well as power in electromagnetic component are not straightforward to estimate. 
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Rawlings and Sauders [55] give estimates of jet power in FR II galaxies based on assumption 
of equipartition between the magnetic field and relativistic electrons and expansion dynam- 
ics of hot spots. According to this minimum energy estimate, powers of FR II galaxies reach 
Xq45 _ iq4:7 c];g/scc, whilc FR I arc an order of magnitude less powerful. This limits the 
possible acceleration sites to FR 11 radio galaxies which are typically far apart. The power in 
strongly beamed sources, like flat spectrum radio quasar, are harder to estimate: radiative 
modehng [56] gives estimates comparable to FR II galaxies, up to 10^^ erg/s. In any case, 
requirements on luminosity (and assumption of high proton content in UHECRs) excludes 
prominent low power nearby sources, like FR I radiogalaxy Cen A, nearest AGN M87 and 
starburst galaxy M82. We favor high power radio galaxies at intermediate distances as 
sources of UHECRs. 

The observed number density of highest power sources is a power law in energy, with 
sources of L ~ lO'^^erg/s having density ~ 10~^Mpc~^. This is smaller by approximately 
an order of magnitude than the density required to produce the small scale clustering, 
~ 10"^ Mpc~^, [72], but, first, the inferred density has a large uncertainty and, second, the 
intrinsic electromagnetic luminosity is definitely higher than the observed radiative one. 

Shape of GZK cut-off. Small spacial density of sources also has a clear prediction for 
the shape of the GZK cut-off [54] . Our model predicts very steep GZK cutoff. 

Total energy budget. The model naturally avoids the energy budget problem posed by 
the steep required injection spectrum (in which case total energy density is determined by a 
somewhat arbitrary lower energy cut-off). In the proposed mechanism, below some limiting 
energies particles just do not leave the jet. Though details depend on jet properties, requiring 
that a cosmic ray crosses a 10 kpc jet over acceleration of jet length of 100 kpc, one can 
estimate Emin ~ 10^^ eV. 

High energy photon and neutrino traces. If high power, fairly distant AGNs are 
sources of UHECR, one may expect that as particle travel towards the observer they initiate 
electromagnetic cascades due to photo-pion production on CMB. Most of the energy of 
electromagnetic cascades ends up in photons in the GeV-TeV range. If propagation is almost 
rectilinear, one may expect an increase of high energy photon flux towards a source at a 
level of ~ 10^^'' cm^^ s^^ [75]. This flux can be (barely) detected by Cherenkov experiments 
hke HESS. Deflection of UHECR in locahzed regions of high magnetic field , associated, for 
example, with supergalactic plane may easily make this fiux to be distributed over large area 
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of the sky. To produce considerable deflection a magnetized sheath of thickness of 1 Mpc 
with magnetic field of the order of 10~^ is required [76]. This may also explain non-detection 
of UHECR towards a particular AGN. Finally, significant neutrino fluxes can be generated if 

protons are accelerated up to 10^^ eV energies [77]. We note however that the induced 10^® 
eV neutrinos arc below the currently advertised threshold for EUSO, OWL and ANITA, 
and most of the potential events will go undetected [78] . 

VIII. DISCUSSION 

In this work we have considered a new mechanism of acceleration of UHECRs to energies 
above the ankle by inductive electric fields in sheared flows. The main points are: 

(i) In a sheared flow protons are at a maximum electric potential energy z/B ■ V x v < 0, 
so that kinetic drift motion along the shear will lead to fast energy gains. If a jet is generated 
by rotating system of black hole-accretion disk with the velocity falling off with cylindrical 
radius, then the quantity B • V x v will be determined by the product B • Q at the source, 
where is the angular velocity of the black hole-accretion disk system driving the jet. 

Radial drift motion may occur if a jet is inhomogeneous along its axis, so gradient drifts 
will lead to energy gains. Alternatively, one may imagine that the toroidal magnetic field is 
perturbed azimuthally and is sheared by differential velocity flow. This will generate non- 
axisymmetric poloidal magnetic field , which would result in radial drift. Finally, particle 
diffusion, resonant or due to field line wondering, may lead to radial particle motion. In any 
particle must cross the magnetic field fines in order to gain energy. 

(ii) Since drift velocity is an increasing function of energy, acceleration rate is also an 
increasing function of energy. Thus, from a pre-accelerated population this mechanism will 
pick up the highest energy particles and boost their energy further. In addition, since energy 
gain is independent of the direction of the drift, pre-accelerated population may be located 
outside of the jet and be puUed-in by gradient drift. Pre- acceleration may be achieved, for 
example, by non-linear shock acceleration in supernova and/or cluster shocks [79]. 

(iii) After a particle has crossed a substantial fraction of the total potential available in 
the system, its Larmor radius becomes comparable to the flow scale. In this case exponential 
acceleration may proceed without gradient drift. As a result, the acceleration rate reaches 
the theoretical maximum (acceleration on time scale of one gyration) and finally a particle 
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can leave the system. 

(iv) Since acceleration is fastest for highest rigidity particles, from a population of a pre- 
accelerated particles with mixed composition, particles with smallest charge - i.e. protons - 
are accelerated most efficiently. We thus predict a dominance of protons over heavy nuclei 
for energies above the ankle. 

(v) Drift acceleration can also produce hard spectra, which are required if sources are 
located at the edge of the GZK sphere. We have also calculated the spectrum of UHECRs 
assuming impulsive injection of pre-accelerated particles with a steep power law spectrum. 
Due to inductive acceleration the energy spectrum flattens and reaches asymptotically a power 
law index of 2. 

The best astrophysical location for operation of the proposed mechanism is cylindrically 
collimated, high power AGN jets. The proposed mechanism cannot work in spherically (or 
conically) expanding outflows since in this case a particle experiences polarization drift, 
which is a first order in Larmor radius, due to the fact that in the flow frame magnetic field 
decreases with time. For a constant velocity flow this drift is always against the electric field 
(for a positively charged particle) and lead to the decrease of energy on time scale R/{cr), 
where i? is a distance from the central source and F is the Lorentz factor of the flow. [A flow 
must acceleration faster than F ~ r in order to beat this polarization drift.] Comparing with 
acceleration time scale (32), and assuming Le ~ Lb ~ -R/F, we conclude that the expansion 
time scale is shorter than the acceleration time scale by a factor ~ ^JRc/^Vujb ^ 1- Thus, 
in spherically expanding flows adiabatic losses always dominate over regular energy gain due 
to drift motion: the proposed mechanism would fail then. On theoretical grounds, AGN 
jets (or at least their cores) may indeed be asymptotically cylindrically collimated [44]. 
Observations of large scale jets, e.g. Pictor A, do show jets that seem to be cylindrically 
collimated on scales of tens of kiloparsec. 

In cylindrically collimated parts of the jet acceleration can happen from sub-parsec to 
hundreds of kiloparsec scales: as long as the motion of the jet is relativistic the total electric 
potential remains approximately constant. Thus, UHECRs need not be accelerated close 
to the central black hole where radiative losses are important [cf. 11]. After a jet has 
propagated parsecs from the central source radiative losses become negligible, see Eq. 6. 

Another constraint on the mechanism comes from the requirement that in order to pro- 
duce radial kinetic drift the jet magnetic field should be inhomogeneous along the axis. 
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Though shocks provide possible inhomogeneity of magnetic field ((^-function inhomogeneity 
on the shock front), we disfavor shock since particles are advected downstream and can- 
not drift large distances along shock surface. In a gradual inhomogeneity a particle drifts 
orthogonally to the field gradients and thus generally will remain in the region of inhomo- 
geneous fields. [One may imagine a hybrid model, where scattering brings a particle back 
to shock front, similar to Eichler's model of shock acceleration near Earth bow shock, [20]]. 
Extragalactic jets are expected to have axial inhomogeneities, both due to non-stationary 
conditions at the source and due to propagation of compression and rarefaction waves gen- 
erated at the jet boundary via interaction with surrounding plasma. 



IX. PREDICTIONS 



Our model has a number of clear predictions, some of which are related to astrophysical 
association of acceleration sites of UHECRs with AGN jets [2] and some are specific to the 
model: (i) one needs a relatively powerful AGN, with luminosity > 10^^ erg/sec. This lim- 
its possible sources to high power sources like PR II radiogalaxies, radio loud quasars and 
high power BL Lacs (fiat spectrum radio quasars). Powerful AGNs are relatively rare and 
far apart, so that a steep GZK cut-off corresponding to large source separation should be 
seen, (ii) UHECRs come from sources with low spacial density. This may be reflected in 
the distribution of arrival directions, (iii) Extragalactic UHECRs should be dominated by 
protons, (iv) Depending on "extra-galactic seeing conditions" arrival directions of UHECRs 
may point to their sources, though complicated magnetic field structure may erase this cor- 
relation. In addition, the fact that only fiows with negative shear can accelerate protons 
implies that only approximately half of such AGNs can be sources of UHECRs (this assumes 
that the AGN central engine - black hole or an accretion disk - is dominated by large scale, 
dipolar- like magnetic field ) . 
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we find j = divE v, which gives A$ = ^(i^^/cV^) • 

APPENDIX A: IMPOSSIBILITY OF ACCELERATION OF UHECRS IN DC 
ELECTRIC FIELDS 

Acceleration by DC electric fields, when there is component of electric field along magnetic 
field or when E > B, is fast and very efficint way of gaining energy. On the other hand, 
in this appendix we argue that UHECRs cannot he accelerated hy DC electric fields. There 
are two reasons. First, if we have a gap (region where there is an uncompensated electric 
field parallel to magnetic field), any lepton will undergo acceleration and after passing a 
potential drop of the order of lO*^ — 10*^ eV will produce a pair, e.g. through inverse Compton 
effect or curvature emission. The secondary pairs will similarly be accelerated, producing 
an electromagnetic cascade that will short out the initial DC electric field . Secondly, even if 
one manages to create an extremely photon-clean surrounding with straight magnetic fields , 
the electric field will accelerate both signs of charges, producing field-ahgned current. The 
magnetic field, associated with these currents would have prohibitively large energy as we 
show below. 

Consider a one dimensional double layer with potential difference <^tot that can freely emit 
particles from each end. To account for possible effects of rotation and associated charge 
density we introduce Goldreich- Julian density naj in the potential equation The governing 
equations are 

9^0 = 4:7r{niei - n^e + noje) 

je — erieVe — COUSt 

ji = CiUiVi = const 




1 m^c' = eii^tot - 4>) (Al) 



where z is the coordinate along the electric field, n and e are density and charges, j is the 
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current density, v is velocity of two species. Thus 

y/(J)j2meC^/e + (j) 

■^e = 2/ I A ' 



= \ c (A2) 

Putting into potential equation, we find 

,// 47r / je(e0 + mec2) j^^^^c^ + ei(<l>toi - 0)) \ 

= — noje , H — , = (,Adj 

c \ ^e0(e0 + 2mec2) ^ei($tot - (/))(2mic2 + ei($tot - 0))/ 

which can be integrated once. The solution of (A3) should satisfy the boundary conditions: 
at 0(z = 0) = 0, 4)'{z = 0) = 0, 4>{z = //) = $tot, 4>'{z = //) = (there is no electric field 
on each end). Thus the problem is an eigenvalue problem on currents. We find 



2mec7e + ^tot ^ ngjejC 

jeCi 2miC^/ei + ^tot ji^2mec'^ /{e^) + 1 

which generalizes the relativistic Child's law [80], see also [27]. Equation for 4>' then becomes 
0'^ 47rje / r — / f— V^tot - (pMrnid^/ei + ^tot - 0' 



I y ^;/i/gi^ / I -^tot 1 V tot I ; 

c V V y/2micyei + ^ 

-^^2m,cye + - Utot - - ^ 



$tot($tot - 0)(2m^cVe^ + 'I'tot - 0) 



2miC^/ ei + %ot 



A5) 



dimensionalizing me = /J^rrii, ^ (f>^tot, ^tot — IqI^iC^ / {Ze) (70 has a meaning of a maximal 
Lorentz factor that ions can reach), noj — > Pajji/iec) and z — > lo^^ yf^^f^, equation for 
0' becomes 



In the limit of highly relativistic ions, 70 ^ 1, this becomes 

0^^ ^ + (Z + /3G^)/i(l-0) ^ ^ 

2 ^/7b ~^/7b ^ ' 

This determine the length of the double layer in dimensionless units. 

z = v^7o^' (A8) 
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Returning to dimensional notations we find that the maximum potential a double layer of 
size R, carrying a current / (so that ji ~ I/2R^, assuming that the width of the current 
layer is of the order of its length) may have is 

^tot = V2^\l^ (A9) 

The maximum energy is then 

~ ^Jlm^cZe (AlO) 
Since / < Etot/ R, where Etot is total energy in a system, 

< (^) (^) (^) (All) 

Thus, the maximum energy achievable by a particle scales only as E^/'^ (c.f. with maximum 
energy achievable in inductive field, oc L}/"^, Eq. (4)). Even the largest energy reservoirs, 
lobes of radio galaxies, cannot provide large enough potential for DC acceleration of UHE- 
CRs. For magnctars with B ~ lO"*^^ G and R ~ 10^ cm estimate (All) gives £max ~ 1.6 x 10"*^^ 
eV, while for supermassive blackholes with B ~ 10"^ G and R ~ 10^^ cm, £max ~ 1-6 x 10^^ 
eV. This excludes the possibility of UHECR acceleration by DC electric fields in astrophysical 
environment. 

There is only one way to circumvent this argument: acceleration region must be really 
devoid of any plasma [as, for example, suggested by 11, 81, in the immediate surrounding 

of dead quasars] . 

APPENDIX B: CAN GRB SHOCKS ACCELERATE UHECR? 

Simple estimates (4) shows that only very luminous sources have enough potential to 
acceleration UHECR. GRBs, with peak luminosity of the order of 10^^ erg/s and total 
potential ^grb ~ 4 x 10^^ V, are then prime suspects [eg 82]. Still, we disfavor this 
possibihty for the following reasons. Internal shocks, where UHECRs are supposed to be 
accelerated, are short transient events occurring at distances ~ 10^^ — 10^^ cm. In order to 
overcome radiative losses, from Eq. (6) it follows that Lorentz factor should be F > 300. The 
first problem is that for spherically expanding ultra-relativistic fiows the available potential 
is smaller than (4) by a factor l/(rA^) ^ 1, since particle can cross sideways only a distance 
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R/r. To see this, consider a conically expanding flow with total luminosity L and opening 
solid angle AQ — ttA^^. [All the relations below apply, roughly, both to large scale and 
small scale magnetic fields.] If b is the magnetic field in the fiow rest-frame, using (3) we 
can write the rest- frame magnetic field as 



a 



V 1 + (7 V c rA0R 
The available potential in the rest frame is then 



(Bl) 



while in the laboratory frame it is 

(B3) 

(since both electric field and magnetic fields are larger in laboratory frame by a factor F). 
At the time when internal shocks are supposed to occur, FA^ 1, so that the available 
potential (B3) is smaller than the total potential by a factor FA^ ~ 30 (for F = 300 and 
A^ ~ 0.1). Potential (B3) is still large enough, so that sub-equipartition fields, a > 0.01 are 
sufficient. 

[In passing we note that the restriction discussed above applies only to UHECR ac- 
celeration at GRB internal shocks occurring at small distances. Since the total potential 
remains approximately constant as long as the expansion is relativistic, F > 1, but the avail- 
able potential increases for decreasing F until F ~ 1/A^ and remains constant after that, 
it is feasible that UHECRs are accelerated at later stages, when the fiow becomes trans- 
relativistic. This cannot, though, be done neither at external shock since magnetization 
there is low, nor at internal shocks, since the fiow becomes trans-relativistic long time after 
internal shocks have been dissipated.] 

Secondly, and most importantly, as long as particle has not crossed all the available 
potential (B3), its Larmor radius is smaller than the size of the fiow in its rest frame 
vl < R/T. After crossing fraction ^ of the total available potential in rest frame (B2) the 
maximum Larmor radius becomes 

rL = i'^ = 4' 
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In a framework of shock acceleration, if Larmor radius is smaller than the size of the system, 
^ < 1, a particle is advected with the flow and is then subject to strong adiabatic losses, losing 
all the energy that it gained from acceleration process. Thus, in order to avoid adiabatic 
losses a particle must cross all the available potential so that its Larmor radius would become 
of the order of the flow scale and adiabatic approximation would break down. Condition (B4) 
is equivalent to the requirement that acceleration time scale for highest energy UHECRs is of 
the order of flow expansion. Thus, the paradigm of shock acceleration of UHECRs in GRBs 
requires that particles are accelerated on time scale of inverse relativistic gyro-frequency, 
nothing less will do. We consider this condition, though formally within the hmits of shock 
acceleration, as a very restrictive and thus disfavor UHECRs at internal shocks in GRBs. 

APPENDIX C: ELECTROMAGNETIC WAVES IN TRANSVERSELY 
SHEARED PLASMA 

In this Appendix we consider properties of electromagnetic waves in transversely sheared 
cold plasma. The results presented here will server as a basis for background electromagnetic 
field that lead to particle acceleration. One of the key differences between our approach and 
previous works [34] is that waves in transversely sheared plasma cannot be considered self- 
consistently in the framework of non-relativistic MHD, since, as we showed in Section HI 
shear introduces charge density, which is neglected by non-relativistic MHD. In addition, 
for perturbations propagating orthogonally to magnetic field , inertial Alfven waves play 
a crucial role, so that induction current, which is neglected in non-relativistic MHD, is 
important. 

Consider cold plasma composed of two species a = e,p. Maxwell equations, continuity 
and equations of motions then read (in this Section we assume c = 1) 

dtB = —curl E 
curlB = 47rJ + at£; 
div = 47r ^ riaCa 
div 5 = 

dtna + V{vna) = 
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dtVa + (Va • V)V« ^^{E+{VaX B)) (Cl) 

rria 

Below we assume that electrons (subscript 1) carry a charge — e while positive chargers are 
protons with charge e. 

Assume that the magnetic field is along x-axis — Bq and that there is sheared fiow 
along y-direction with V — Vz{y). Then, there is electric field Ey — —Vz{y)Bo and corre- 
sponding charge density 

Po = ^divE=-^V'Bo (C2) 

Consider small fiuctuations on this background oc ex.p{i{kxX + kzZ — ujt)}, B = Bq + 56, 
E = Eq + 5e, ria — rio + ria, v = V + v^. Assuming, for simplicity — 0, perturbations are 
governed by equations 

rj kz^z r, ^6 ie + kz€y 
dby = , dbz = — , dba: = 

UJ 00 UJ 

iuty + A'Knoe{viy - V2y) = 

iujCz - Sb'^ = 47rnoe((ni - n2)V + no{viz - V2z)) 

5b'^ — i{kz5by — cucx) 

5b'y = -ikzSbz 

e'y + tkzCz = 47re(ni - 712) 

—LOUa + kzUoVix — in^v'^y + kzUiV = 

i^Ba^z + i iy^DVaz + Vay{uJBa + V')) =0 

= 5byV (C3) 
where a;^) = a; — /c^V. Eliminating 5b and Ua- 

= ^^^"^ ~ '^"^ (C4) 

we derive the velocity fluctuations as 

Vaz n n , , ^ ^ 

where Da = ujIq — uisa — '^BaV- Equation for Cx separates 

4 = {u^' - kl)ex, (C6) 
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while equation Cy and Cz obey a system of equations 

= ~2 U (4vrenoa;(viy - V2y) + fc^e'J 

LO — 

ik^e^ + 4 ^ ° (/^^('yu - ^^2z) - - ^;22/)) (C7) 

where velocities can be eliminated using (C5). In full generality the system becomes too 
complicated to be reproduced here. We consider two simple case: propagation along y-axis 
and a low frequency limit cu <^ cusa- 

1. Propagation along z axis 

For wave with = 0, we can immediately find dispersion relation for the fast mode 
(inertial Alfvcn mode) polarized along y: 

u^^kl + '-^i- ^%^] (C8) 

We are interested in low frequency waves ui <S ujbi, ojb2- Assuming that plasma is 
electron-ion, (first type of particles are electrons with |a;si| 3> a;s2), and introducing mag- 
netization parameter 



we find 



a = (C9) 



a;2 = kl + '^ (CIO) 



In the absence of shear this gives dispersion relation for inertial Alfven waves uj — 
'JaJ(T+~a)kzC — VAkz where Va — y^ct/I -|- ac and we have reinstated explicitly the speed 
of hght. This relation is relativistically exact; in the non-relativistic hmit cr <S 1, this gives 
so called inertial Alfven waves with Va — \foc — B j ^JA.'Km^nQ. 

In the non-relativistic hmit Va, 1^ <S 1 Drift dispersion relation (CIO) gives 

u^K{yA±v) (cii) 

2. Oblique propagation normal to the magnetic field : low-frequency limit 

Next we consider low frequency waves propagating obliquely to velocity shear. In the 

limit ijj <^ujb2 ^ f^si we find 

_ _ ( u:T){kzi^D - (c^^ - A;.)^V)e; - cu(cu^ - A;,^)e,r \ 
V BouUB2{u^^-kl + u:lla) ) 
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Viz - V2z = ey = — — — (C12) 

Bou}B2 o- uj^ - kl + ujjj/a 

which gives an equation for continuous spectrum normal modes 



V<7(a;2 _ ^2) +u;lj ^ {a{u;' - kl) + u^lf 

I A + ^ I g K{K-uV)V'^ \ {K-uV)V" 

Assuming weak non-relativistic shear V (neglecting quadratic terms in V) , this reduces 
to 

f , / l^^^ , 2 + + ^-T^K^r^-^ < = (CM) 

\(t{uj^ - kl) + ujI ' j \ a ujdHuj^ - kl) + ujI) J 

This equation is reminiscent of Rayleigh equation in compressible fluid [83] : second derivative 
of the velocity profile is divided by the shifted drift frequency. This determines the usual 
Rayleigh instability. (For a flow to be Rayleigh unstable, the necessary condition is presence 
of inflection point where V" = 0.). [Also note, that for kz = Eq. (C14) gives uj = VAky, so 
that shear does not affect wave propagation in this case.] 

For linear profile V" — and in the non-relativistic limit V <^ 1, Eq. (C14) reduces to 

Yl 

~ k^ V 

Far from the resonance points equation of the form 



8, (mo + = 

with slow varying function /(|/) can be solved by WKB method, which gives 

ez oc --^exp{±i j dy/^Jf{y)} (C16) 

Properties of small perturbations on the background flow are closely related to the sta- 
bility of the flow. Ideal sheared flows may be subject to Kelvin-Helmholtz and Rayleigh 
instabilities [83]. In addition, in warm plasma new types of instabilities appear, which are 
related to ion cyclotron motion [34]. Discussion of the stability properties of the transversely 
sheared flows is beyond the scope of the paper: we assume that the sheared configuration 
under consideration is stable. This can be shortly justified by noting that Rayleigh-type 
instability depends on the presence of inflection point in the flow and can be eliminated by 
choosing, e.g. , linear velocity profile, while kinetic instabilities are negligible if the back- 
ground plasma is cold. 
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FIG. 1: Potential energy of a charge in a flow for positive and negative shear. Cylindrically 
collimated flow carries toroidal magnetic fleld5<^ and is sheared, Vz{r), so that regions closer to 
the axis move with larger velocity. There is an inductive electric field, directed either away from 
the axis (for positive shear) and towards the axis (negative shear). In a rest frame moving with 
the flow at some r = r(ro) the electric field is either towards the origin (for positive shear) or away 
from the the origin (negative shear), so that a proton is either at potential minimum or maximum. 
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FIG. 3: Particle trajectory y(t) in strongly sheared layers with Vz{y) = r]y and for different values of 
■q measured in term of the Larmor frequency lob = eB^/{mc). Magnetic field = Bqx is constant. 
A particle initially at y = 0, with velocity along y axes and Lorentz factor 70 = 3. For 77 < ricru = 
y^7o — 1 — 7o = —0.171 trajectories are unbound and particle energy increases approximately 
exponentially, while for larger ij trajectories are bound. Particle displacement is measured in term 
of non-relativistic Larmor radius c/ojb- For rj = rjcrit particle reaches ycrit = '^/'^B\J^/Vcrit ~ 
moves along z axis with constant Lorentz factor jcrit = 7o(l + Verity crit^B / {'^(^)) ■ 
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FIG. 5: (a) Maximal energy gained by a particle in a non-linear inertial Alfven wave propagating 
in a sheared flow as a function of initial phase of a test particle for the following parameters: 
Va = 0.5c, a = 0.25, kz = 0.025, Us = —5 (both in units c/ojb), initial energy 70 = \/2. Upper 

curve corresponds to the full case shear plus a wave, lower curve is the case without shear, (b) 
Dependence of maximal energy gain on shear scale, ys = —3, ys = —10, = — 15 (top to bottom). 
In smaller shear (larger negative ys) acceleration is less efficient. 
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